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We perform renormalization group transformations to construct optimally local per- 
fect lattice actions for free scalar fields of any mass. Their couplings decay exponentially. 
The spectrum is identical to the continuum spectrum, while thermodynamic quantities 
have tiny lattice artifacts. To make such actions applicable in simulations, we truncate 
the couplings to a unit hypercube and observe that spectrum and thermodynamics are 
still drastically improved compared to the standard lattice action. We show how precon- 
ditioning techniques can be applied successfully to this type of action. We also consider a 
number of variants of the perfect lattice action, such as the use of an anisotropic or trian- 
gular lattice, and modifications of the renormalization group transformations motivated 
by wavelets. Along the way we illuminate the consistent treatment of gauge fields, and we 
find a new fermionic fixed point action with attractive properties. 
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1 Introduction 



Renormalization group transformations (RGTs) are an important operation, both, in field 
theory and in statistical mechanics. They allow to vary the cutoff of a system without 
altering its physical contents. In particular RGT fixed points and their vicinity are of 
interest, both, conceptually and for practical purposes. 

If a system is defined on a (Euclidean) lattice, then RGTs can be performed by block 
variable transformations and a subsequent rescaling |j^] . In particular for the simple system 
of a free scalar field, it has been shown a long time ago that such block variable RGTs - for 
instance with blocking factor 2 ~ can be iterated an infinite number of times in momentum 
space 0. For suitably chosen RGT parameters and massless scalars, one arrives at a finite 
fixed point action (FPA). There the observables are independent of the lattice spacing - 
which is an inverse UV cutoff - so that we can represent exact continuum physics on the 
lattice. 

In each transformation step the average of the variables in one block of the fine lattice 
is related to the corresponding block variable on the coarse lattice. This has often been 
implemented hy a 6 function. However, the 6 function is just one possibility in a class of 
implementations, which keep the partition function and all expectation values invariant. 
In Ref. a generalization was suggested, where this relation is smoothly implemented by 
a Gaussian, and it was observed that in this way the locality of the FPA can be improved. 
Q This is very important for practical purposes, because one ultimately has to truncate the 
couplings to short distances, and to hope that this truncation does not affect the perfect 
properties too much. If one sends the blocking factor to infinity, then one single RGT 
leads to the FPA |3|, This amounts to a technique that we call "blocking from the 
continuum" : one starts from a continuum theory and defines lattice variables by relating 
each of them to the Riemann integral over a lattice cell. 

A similar construction has been realized for free fermions with a 6 function RGT 
and a generalized RGT Q. In that context the FPA also yielded new insight into the 
fermion doubling problem. By a somewhat more complicated RGT one obtains a FPA for 
staggered fermions too Also here we can optimize locality by using a smooth RGT |p 
and the result can be reproduced again by blocking from the continuum |^, IC]. Moreover, 
blocking from the continuum also leads to a fixed point for free non-compact gauge fields 
0- 

FPAs have the remarkable property that a system is regularised but not contaminated 
by any cutoff artifacts. Of course it is a dream to construct such a regularisation also 
for interacting theories at arbitrary correlation length ^ (FPAs only exist at ^ = oo). 
Indeed, general arguments show that such "perfect lattice actions" do exist even at finite 
correlation length Q, but they are hard to construct. By blocking from the continuum this 
can be achieved at least perturbatively, and the perfect action was constructed explicitly 
to first order in the fermion-gauge coupling for the Schwinger model and for QCD 

^Also the scalar FPA of the S function RGT is local in the sense of an exponential decay of the couplings 
(in contrast to the fermionic FPA), but for certain Gaussian RGTs this decay becomes even faster, see 
Section 2. 
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in], P, pi. For the anhar monic oscillator, the first order perturbatively perfect action 



has been constructed and simulated in order to test its scaling and asymptotic scaling 
behavior Q. 

This approach to fight artifacts due to the finite lattice spacing - which are the worst 
systematic errors in many Monte Carlo simulations - has picked up new momentum since 
the appearance of Ref. ||l5|. There the authors noticed that for asymptotically free theories 
the construction of a FPA (with rescaling of the weakly relevant coupling) is a classical 
field theory problem. Q Hence it can be solved by minimization instead of (numerical) 
functional integrals, which simplifies the task enormously. Moreover the authors suggest 
to use the FPA, which is a classically perfect action, even at moderate correlation length, 
and they demonstrated for the 2d 0(3) model that this can be a very successful approxi- 
mation of a perfect lattice action. The construction can be performed non-perturbatively, 
by numerical inverse blocking RGTs. Of course for complicated theories like QCD, or 
pure SU{3) gauge theory [0], even this approximation is tedious. In addition, problems 
related to the parameterization and truncation of such quasi-perfect actions are difficult 
to handle, but crucial for potential applications. With this respect it is advisable to start 
from a simple situation and study there, which procedures are promising. The present 
paper is a contribution to that issue. A synopsis of our results was anticipated in Ref. |18]. 



In this paper we deal with scalar fields. In analogy to the previously studied fermions, 
the free particles provide an improved formulation, which is promising also in the presence 
of interactions. In Section 2 we derive perfect lattice actions for free lattice fields by 
blocking from the continuum and we optimize their locality. In Section 3 we truncate the 
couplings to a unit hypercube and consider the effect on the dispersion relation. In Section 
4 we illustrate the improved thermodynamic scaling properties of this "hypercube scalar". 
For the evaluation of the scalar matrix - analogous to the fermion matrix - we can still use 
preconditioning in the hypercubic case, as we show in Section 5. In Section 6 we discuss 
perfect actions on anisotropic and triangular lattices. As a further alternative we consider 
in Section 7 a class of new blocking prescriptions, which is motivated by the wavelet resp. 
B-spline formalism. The most promising new variant is also applied to fermions (Appendix 
B). Section 8 contains our conclusions and an outlook on applications. 



2 Perfect lattice scalar actions 

Let us start from a scalar field <j) and its action S[(j)] on a hypercubic lattice in d dimensional 



Euclidean space. The Kadanoff transformation to a coarser lattice can be written as |19| 



^Alternatively there are many attempts to use real space blocking RGTs (without classical approxima- 
tion). For recent work on scalars resp. pure SU{3) gauge fields, see Ref. [|l6|. A limitation there is that 
one can hardly implement the blocking constraint in a way different from the 5 function. 
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where S'[4>'] is the action of the new scalar field (p' on the coarse lattice. We integrate over 
the lattice measure T>(j) = jY\xd(t)x- The kernel K[(j)'^(j)\ has to be chosen such that the 
partition function and its derivatives remain invariant under this transformation (up to a 
possible constant factor). This can be achieved by the requirement 

ypc/.' K [(/.',(/>]= const. (2.2) 

The simplest and most popular choice is 

K[(P',4>]=\{5{cP',,-Ax,[ct^]). (2.3) 

x' 

Here x' runs over the coarse lattice, and the symbol A^i [cp] is an average taken in the block 
of the fine lattice, which is attached to x' . For instance, one may just take the arithmetic 
mean value together with a rescaling factor, 

^-'M = ^E'^- (2-4) 

This sum extends over all the n"^ fine lattice points in the block associated with x' . After 
the transformation we express all quantities in units of the coarse lattice. The factor 6„ is 
used to neutralize this rescaling so that we can arrive at a finite FPA. Hence its value is 
determined by the dimension of the scalar field, 

6„ = n('^-2)/2. (2.5) 

However, the choice ( |2.3| ) is by no means unique. A class of generalizations has been 
considered by Bell and Wilson, where the 5 function is replaced by a Gaussian [^, 

i^[</<',0] =exp{--5][</>^,-^,,</.]2} , (a>0). (2.6) 



In the limit a ^ we return to the 5 RGT ( p. 31) , but the condition (2^) holds for all 
positive a. 

In Ref. Q the fixed point for a massless scalar particle has been computed by an 
infinite number of iterations of such block factor 2 RGTs. The iteration starts from the 
standard lattice action in momentum space, 

'^['^^ = (2^/"^'^ c/>(-A;)^p0(A;) , (k, := 2sin^), (2.7) 

where we set the lattice spacing equal to 1. 

However, instead of this tedious iteration we can also send the blocking factor n — ^ oo, 
so we obtain a FPA by performing only one RGT ^. Since we call the spacing of 
the final lattice 1 again, the original lattice variables move closer and closer together (in 
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the final units) as n increases. Finally the sum in condition ( p. 2D turns into a Riemann 
integral, 



A,ip= / Avp(y). (2.8) 

J Cx 

Here we do not start from a field (px on a fine lattice any longer, but from a continuum 
field ^p{y)- We integrate over hypercubic unit cells Cx with centers x, and relate those 
integrals to the new lattice variables (px- Also the action on the fine lattice is replaced by 
the continuum action. Now the RGT step requires a continuum path integral, so it can 
only be computed in certain cases. Of course this limit n — > cx) can also be taken for the 
generalized form (12? 



We first reproduce this calculation in its short-cut form, including also an arbitrary 
mass m of the scalar field. There is no FPA at m > 0, but we construct an entire 
renormalized trajectory, i.e. a curve of perfect lattice actions, for free scalars of any mass. 
It is parameterized by m and crosses the critical surface in a fixed point at m = 0. 

The lattice action is given by 



(2-kY Jb 

k - 2TTl)[{k + 27Tlf + m'^](f{k + 2t:1) 



+i(j{-k)[^{k) - ^ ^{k + 2TTl)Ii{k + 2-kI)] + %{-k)a{k)\ }, (2.9) 



where B = ] — vr, vr] is the (first) Brillouin zone. We have introduced an auxiliary lattice 
scalar variable a, and the function 11 is defined as 

ii{k) = n 1^. (2.10) 



We can write eq. (2^) as a convolution. 



A^^^j.l'yn.-yMy). /(„) = { J I""' "^1= . (2.11) 

and n(A;) is the Fourier transform of f{u). Instead of the piece-wise constant function 
used here one may also insert other functions /(n), if they fall off sufficiently fast. Some 
options will be discussed in Section 7. 

We denote the continuum propagator as 

substitute 

if{k + 27rl) = ip{k + 2-nl) - ia{k)A{k + 27rl)U{k + 2-irl), (2.13) 



5 



and carry out the Gaussian integration over (p: 

e-^[<^] = lvaexp{jl^ j ^d''k[ia{-k)m -\<-k)G{k)a{k)]} , 
G{k) := A(A; + 27r0n2(/fc + 27r0 + a. (2.14) 

By further substituting 

a{k) = a{k) - iG{k)-^ ct){k) (2.15) 

and integrating a we arrive at 

^^^^ ^J^^Ib "^'^ \^^-^)G{k)-^m. (2.16) 

which shows that G{k) is the perfect free lattice propagator. Note that this scalar perfect 
lattice action is also the renormalized trajectory of the 0{N) non- linear a model in the 
limit — > oo (up to a factor N) ||2^, in analogy to the large limit of the Gross 
Neveu model The sum over / G Z'^ converges in any dimension thanks to the 11 

function, but it can be computed analytically only in d = 1. If we now specify the RGT 
"smearing parameter" (because it smears the 5 function to a Gaussian) to be 

sinh m — m 

a = aim) := ^ , (2-17) 

then the perfect action in d = 1 becomes ultralocal, 

sinhm-m^ 1 - n ■ i e^ 

Gd=i{k) = 3 J- — , m:=2smh— . (2.18) 

m-^ k^ +m? 2 

In this case, the couplings are even restricted to nearest neighbors. In particular the Id 
fixed point propagator (defined at m = 0, where a = 1/6) is identical to the standard 
lattice lattice propagator l/k^. This is obviously the choice for a, which optimizes the 
locality in d = 1. For a general blocking factor n, the Id perfect propagator turns ultralocal 
with 

an{m) := a{m) (l - (2-19) 

In higher dimensions - that is, in field theory - locality of the perfect propagator can 
only be achieved in the sense of an exponential decay. For fermions the RGT parameters 
which provide ultralocality in d = 1 have been used successfully also in higher dimensions. 



It turned that they still yield extremely local perfect actions in d > 1 |§, 11, IC]. 
For free gauge fields, a similar (but somewhat more complicated) smearing term leads 
to ultralocality (the standard plaquette action) in d = 2 and extreme locality in d = 4 
|ll| , p^ . Hence we are guided to apply the same strategy to scalar fields and use a = a 
as given in eq. ( |2.17| ) in any dimension. 

Bell and Wilson compared numerically various smearing parameters in d = 3, and 
in the fixed point they found optimal locality at 1/02 — 8, which agrees with the value 
resulting from our analytical method. 
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The perfectness of the action ( p. 16 ) can be verified by an exphcit additional block 



factor n transformation (n < oo). This is demonstrated in Appendix A. The result is that 
a block factor n RGT applied on G{k) at mass m reproduces G{k), now at mass n ■ m. 
This is exactly the expected behavior on the renormalized trajectory. This calculation 
also justifies eq. ( 2.ig| ) as a generalization of a = Ooo given in eq. ( |2.17| ). 



Let us consider the spectrum of the perfect lattice action ( p.l6|) . We denote k = {k, k^i) 
and identify the energy as, E = i{k(i + 2iTlfi). Thus the spectrum is given by 

E"^ = {k + 2TTlf + . (2.20) 

There is a branch for each /. The branch for Z = is the exact continuum spectrum. In 
this sense, the full, continuous Poincare invariance is present in the perfect lattice action. 
We emphasize that the form of the lattice action itself does not reveal this property: for 
instance, the hypercubic structure of the lattice is visible in it, and perfect actions on non 
hypercubic lattices have different forms, see Section 6 for examples. However, physical 
observables do have the full continuum symmetries, without any artifacts due to the finite 
lattice spacing. 

The higher branches in the spectrum are required for perfection, since the lattice im- 
poses 2tt periodicity on the momenta. 

We now address the question of locality. For massless Wilson-like fermions and a 5 
function RGT one obtains a nonlocal fixed point. Its couplings do not decay exponentially, 
but only oc r^~°' |^]. This is the way a contradiction with the Nielsen-Ninomiya theorem 
is avoided. | However, there is no theorem forbidding the scalar 5 function FPA to be 
local. We analyze this issue in d = 1. 



A perfect lattice fermion propagator - constructed in analogy to eq. ( |2.9D - reads |11] 



where is a smearing parameter analogous to a. For aj = and we obtain in one 
dimension 



G^=i(A;)U,=o = - - ^ ctgh^ - ictg^l ' , G ] - 7r,7r]. (2.22) 



This propagator has a zero - indicating non-locality - at /c = vr, iff m = 0. 
The corresponding Id perfect scalar propagator can be written as 

Gd=i{k)\a=o = ^[Gii(A:)U,=o + Gi,(-A:)U^=o]. (2.23) 

^For any a/ > the fermionic FPA turns local Q, but then the chiral symmetry of the lattice action 
is explicitly broken, avoiding again a contradiction with the Nielsen-Ninomiya No-Go theorem. However, 
this only means that the chiral symmetry is not manifest in the lattice action. In physical observables it 
is still present similar to the Poincare invariance, which is present in the spectrum, although not 

manifest in the action. 
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Here we don't find any zeros; in particular 



J^oG'd=i(vr)|a=o 



1 

12' 



(2.24) 



Tliis different behavior of the 6 function FPA for fermions and scalars persists in higher 
dimensions. 

In coordinate space we write the perfect scalar action as 

S[(I>] = 1 E Mx-y)(l>y, (2.25) 

where p{r) is the inverse propagator, i.e. the (inverse) Fourier transform of G{k)~^. Fig. 
H shows how the exponential decay coefficient ci on an axis depends on the smearing 
parameter a, p{i, 0, 0, 0) oc exp(— ci(a)i), at m = 0, m = 2 and m = 4. R 




Figure 1: The decay of the couplings \p{i, 0, 0, 0)| oc exp{— ci(a)i} of the perfect scalars for 
masses 0, 2 and 4- The maxima are very close to the values of a{m), 0.167, 0.203 resp. 



0.364, which correspond to eg. (2.17). 

The figure confirms that a given in eq. ( 2.17| ) is an excellent candidate for optimal 
locality. It also affirms that the maximum moves to a larger a as the mass increases. We 
stay with the neat analytical criterion for the selection of a presented above, and focus 
on the mass dependent value of a given in eq. ( p.l7| ). Fig. |2| illustrates the exponential 
decay at various masses. The decay becomes even faster if the mass increases, as it was 
observed for fermions before This is illustrated in Fig. |3|, where we consider the 

exponential decay \p{i, i, oc exp(— C4(m)i). Especially in the interval m = . . . 2, the 
decay coefficients follow very closely a parabola, 



C4(m) ~ 6.563 + 0.247 



(2.26) 



a- 



The largest couplings in 4d coordinate space for m = 0, 1, 2 and 4 are given in Table 



''The curves in Fig. |l| are Bezier interpolations, i.e. smooth curves which do not pass through every 
single data point. Since the exact curve is not that smooth, the Bezier interpolation is appropriate to 
identify the region of optimal locality. 
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r 


m = 


m = I 


m = 2 


m = 4 


(0000) 


4.54276730 


4.46018784 


4.11732967 


2.58349457 


(1000) 


-0.25656026 


-0.21033619 


-0.12051159 


-0.01815118 


(1100) 


-0.06722788 


-0.05285547 


-0.02707061 


-0.00294191 


(1110) 


-0.02211001 


-0.01677521 


-0.00778413 


-0.00061988 


(nil) 


-0.00828931 


-0.00611485 


-0.00261666 


-0.00015765 


(2000) 


-0.00079938 


-0.00042262 


-0.00002726 


0.00001090 


(2100) 


-0.00058401 


-0.0003G971 


-0.00010471 


-0.00000201 


(2110) 


-0.00027747 


-0.00018443 


-0.00006107 


-0.00000192 


(2111) 


-0.00009663 


-0.00006885 


-0.00002660 


-0.00000098 


(2200) 


0.00022983 


0.00016003 


0.00005719 


0.00000171 


(2210) 


0.00013280 


0.00008745 


0.00002724 


0.00000058 


(2211) 


0.00008482 


0.00005317 


0.00001456 


0.00000022 


(2220) 


0.00005018 


0.00003084 


0.00000804 


0.00000010 


(2221) 


0.00003271 


0.00001947 


0.00000464 


0.00000005 


(2222) 


0.00001011 


0.00000573 


0.00000120 


0.00000001 


(3000) 


-0.00005999 


-0.00003548 


-0.00000842 


-0.00000009 


(3100) 


-0.00002060 


-0.00001177 


-0.00000252 


-0.00000002 



Tabic 1: The largest couplings p{r) of the perfect scalar action in coordinate space for 
masses m = 0, 1, 2 and 4- The table includes all couplings with absolute values > 10~^. 
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Figure 2: The decay of the couplings \p{i,i,i,i)\ in the optimally local perfect scalars for 
masses 0, 1, 2 and 4- The lines are least square fits to \p{i,i,i,i)\ cx exp{— C4(m)i}, with 
C4(0) = 6.515, C4(l) = 6.775, 04(2) = 7.493 and 04(4) = 9.769. 

3 Truncation by periodic boundary conditions 

So far we have worked in an infinite volume, now we consider a hypercubic finite volume 
with periodic boundary conditions. Of course in a finite volume there are no fixed points 
in the proper sense, but there are still perfect actions, similar to the perfect actions at 
finite mass. 

Let us fix a volume V and m = 0. Then we start from a much larger volume, which 
shrinks under a number of RGTs (block factor < 00) just to the fixed final volume V. If 
we now expand the initial volume more and more, such that the number of RGTs gets 
larger and larger, then the couplings in V may converge. Their limit represents the perfect 
action in the finite volume. (Analogously we can fix a finite ultimate mass mV = oo and 
start from a smaller and smaller initial mass, to obtain after a number of RGTs couplings 
that converge to those identified for the final mass in Section 2.) 

The transition from infinite to finite volume is technically easy: we just replace the 
integrals over the Brillouin zone in all the formulae of Section 2 by discrete sums, and the 
rest remains unaltered, including the smearing parameter for optimal locality. We may 
call this a "decimation in momentum space". Of course, the couplings p{r) will change a 
little compared to y = 00; for instance an exponential decay is not possible any more. At 
the boundary of V they have to vanish smoothly, which suggests the use of this mechanism 
as a truncation scheme: we impose periodic boundary conditions in some volume, set the 
couplings to points outside this volume to zero, and finally use this set of couplings in any 
volume. 

This scheme has a number of virtues: the normalization conditions for J2r Pi''') 
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12 3 4 

m 



Figure 3: The decay of the couplings of the perfect scalars - measured by C4 in \p{i, i, i,i)\ oc 
exp{— C4i} - as a function of the mass m. The optimally local perfect scalars follow 
approximately a parabola, and its decay is clearly faster the one obtained from a 6 function 
RGT. 



J2r r'^pir) are precisely fulfilled (in particular for m = these sums are resp. —2c?), and 
the transitions to lower dimensions by summing over the lattice sites in the supplementary 
directions is also exact. Hence the selection criterion for the RGT parameters, which refers 
to the mapping on d = 1, is still sensible for the truncated system. 

All these properties are not obeyed if we truncate by just chopping off the couplings 
in infinite coordinate space outside a given volume. If we start in this way, we have 
to re-adjust the normalizations a posteriori, which is quite ambiguous. Moreover, the 
periodic truncation has a huge practical advantage: certain quantities in perfect lattice 
perturbation theory take a rather complicated form, such as the perfect quark-gluon vertex 
in QCD |11, ll^. It is a difficult numerical task to evaluate them and to transform them 



to coordinate space, where they are ultimately needed. If only a few discrete momenta 
are involved, this task simplifies tremendously. 



3.5 



3 



2.5 




Figure 4: The mass dependence of the hypercubic couplings po = p(0, 0,0, 0) (left), as well 
as pi = p(l, 0,0,0), p2 = p{l, 1,0,0), p3 = p(l, 1, 1,0) and p4 = p(l, 1, 1, 1) (right). 
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We apply the preferred procedure by imposing periodic boundary conditions over a 
distance of 3 lattice spacings. Thus the remaining couplings are restricted to a unit 
hypercube. There are formally 3*^ of them, but they come in just d+1 equivalence classes. 
Furthermore the mapping condition on the ultralocal Id case and J2r f'^p{f) impose three 
constraints on the equivalence classes, such that e.g. in d = 4 there remain only two 
degrees of freedom. We give the full set of d + 1 couplings in d = 2, 3 and 4 for masses 
m = 0, 1, 2 and 4 in Table |^, and we illustrate the smooth mass dependence of the 
truncated 4d couplings in Fig. |^. This truncation range is of practical interest, since it 



has been shown that hypercube particles are tractable in 4d simulations |T3| , pq , 26]. 

One may argue that the really perfect actions of Section 2 (in d > 1) are academic, 
since we always need some truncation for practical purposes. As a first check on how 
much the perfect properties suffer from the 3-periodic truncation, we show the dispersion 
relation of the "hypercube scalar" for m = and m = 2 in Fig. ^, where they are compared 
to the perfect or continuum spectrum and to the spectrum obtained from the standard 
lattice action. As an example, we consider the direction (110) and we observe a dramatic 
improvement, all the way up to the edge of the Brillouin zone. The situation is similar in 
the (100) and the (111) direction. This is a success of the optimized locality in Section 2. 
Note that the truncated perfect energies start with a slight overshoot at small momenta. 
Following Symanzik's improvement program - which is fundamentally different from the 
program discussed here - it would be possible to correct the leading order artifacts oc 
Then the dispersion becomes even better for small momenta, but with respect to larger 
momenta this is not necessarily an advantage. At the edge of the Brillouin zone the curve 
has to bend down due to periodicity, so it can be favorable to start with a slight overshoot, 
in order to follow the perfect curve over a wide range. After all, if one wants to simulate on 
very coarse lattices, then not only the momenta |/c| <C 1 are important. We also emphasize 
that our "hypercube scalar" follows directly from a very general prescription; no particular 
tuning for a good dispersion is involved. This raises hope that the same prescription is 
also successful for other quantities. 

For comparison, we consider a Symanzik improved scalar by including additional cou- 
plings along the axes over two lattice spacings. The action 



is O(o^) improved; the remaining artifacts for the free scalars are of O(a^). The dispersion 
relation of this Symanzik improved scalar is also shown in Fig. ^ (left). By construction 
it is close to the continuum scalar for small momenta, but at ki = k2 = 1.239 it is hit by 
a higher branch, and the continuation of the curve is just the real part of two complex 
conjugate solutions, which are useless. 

In exactly the same way one can improve the staggered fermions and arrives at the 
Naik fermion |2^, which combines the discrete derivatives to nearest neighbors and over 
distances of three lattice spacings with the relative weights 9/8 and —1/24. A similar 
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construction starting from the Wilson fermion is called D234 action |28]. It turns out that 
its behavior is very similar to the one of the Symanzik improved scalar shown in Fig. ^ 
(left). Also Naik fermions and D234 fermions are hit by higher branches at \k\ = 0(1), 
which marks the end of the real solution for E{k). Those dispersions are compared to the 
truncated perfect fermions in Ref. 



10, 1^ 



Similarly, if we truncate an optimally local, perfect free gauge field to a 4d hypercube, 
we find an excellent (transverse) dispersion relation jl^ 
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Figure 5: The dispersion relation for a scalar of mass m = (left) and m = 2 (right) in 
the (110) direction with the perfect, truncated perfect ( "hypercubic" ) and standard lattice 
action. For the hypercube action (as well as the standard action) there are no "ghosts" 
(higher branches) in any direction, in contrast to the Symanzik improved action. 



4 Thermodynamic properties 



In thermodynamics the artifacts due finite lattice spacing are particularly bad |29|, hence 
the use of improved actions is strongly motivated. Thermodynamics with quasi-perfect 
lattice actions has also been studied for pure SU (3) gauge theory [30|, and for the 2d 0(3) 
model 1 31]. 

First we consider massless scalars with Nt lattice points in the Euclidean time direction. 
We impose periodic boundary conditions in the 4 direction over those Nt lattice spacings. 
The pressure p on an infinite lattice with the standard action is given by |32] 



p_ 

rp4 



(2^ 
1 

~2^ 



Nt 



d^k 



— lnA(fc, A;4,„) 



lA:4_„=27rn/A''t 



n=l 



dki lnA(A;,A;4) 



(4.1) 



where T is the temperature. 

To study the improvement of the lattice actions from Section 2 and 3, we just replace 
the standard lattice propagator in eq. ([4.1|) by the the fixed point or truncated fixed point 
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(000) 


4.02729212 


4.03823126 


3.87615269 


2.54721313 


(100) 


-0.39376711 


-0.31766913 


-0.17499188 


-0.02402874 


(110) 


-0.11305592 


-0.08745023 


-0.04296647 


-0.00419173 


(111) 


-0.03850230 


-0.02898696 


-0.01310490 


-0.00094393 


d = 4 










(0000) 


4.54224606 


4.45988071 


4.11725747 


2.58349386 


(1000) 


-0.25747697 


-0.21082472 


-0.12055239 


-0.01814036 


(1100) 


-0.06814507 


-0.05342220 


-0.02721974 


-0.00294419 


(1110) 


-0.02245542 


-0.01701401 


-0.00787336 


-0.00062377 


(1111) 


-0.00802344 


-0.00598647 


-0.00261577 


-0.00016008 



Table 2: The couplings p{r) of the truncated perfect scalar action in coordinate space for 
masses m = 0, 1, 2 and 4- 

propagator. In Fig. ^ the results for the three lattice actions are compared to the Stefan 
Boltzmann law in the continuum, 

^ - - (4 2) 

In this respect even the "perfect actions" contain lattice artifacts. The reason is that 
we ignored "constant factors" when we performed the Gaussian functional integrals over 
Dx and Da in Section 2. For instance in the spectrum such factors have indeed no influ- 
ence. However, they can depend on the temperature, yielding artifacts in thermodynamics. 
These artifacts are exponentially suppressed, and - as the Figures show - they disappear 
very fast as Nt increases. The truncation slows down a little the convergence to the Stefan 
Boltzmann value, but the behavior is still strongly improved compared to the standard 
lattice action. The latter suffers from quadratic artifacts (asymptotically 0.38it'^ /N^). Hy- 
percubic truncation also yields quadratic artifacts, but only of about —O.Obn'^/N^. 

As a second example, we return to T = but introduce a chemical potential /i according 



to the instruction for general lattice actions given in Ref. |33|. We then measure the scaling 
ratio p/fi'^ for massless scalars. In the continuum it amounts to l/487r^, and all lattice 
actions approach this value in the limit ^ 0. For finite fx the lattice artifacts are visible. 
Fig. ^ (right) shows that they are dramatic for the standard action, whereas the artifacts 
in the hypercube action are again quite modest. 
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Figure 6: Left: The ratio p/T'^ at Nt discrete points in the Euclidean time direction for the 
perfect, the truncated perfect and the standard lattice action. Right: the scaling ratio p/fJ.^ 
(fi: chemical potential) atT = Oin the continuum, compared to different lattice actions. 

5 Preconditioning the hypercube scalar 

If we want to simulate a hypercube scalar action, we are interested in efficient ways to 
compute tlie scalar determinant, analogous to the fermion determinant. There are a 
number of well established algorithms for this purpose, for an overview see Ref. |34]. 
The maximal speed of such methods is limited by the eigenvalue distribution. The goal 
of "preconditioning" is to transform the matrix such that its spectrum becomes more 
favorable with this respect. In particular, it is of advantage to squeeze the eigenvalues 
into a narrow interval. 

In the case of the standard action we can order the matrix lines by passing successively 
through the even and the odd sub-lattice, and we obtain a matrix of the block form |^ 



M 



II -A 
-A 1 



(5.1) 



In a volume V , M is aV xV matrix which splits into 4 blocks of the same size (1 represents 
unity). We decompose M into M = 1 — L — [/, where L, U are the strictly lower resp. 
upper triangular part of M. Now we define the matrices Vi := 1 — L, V2 ■= Tl — U , which 
can be inverted trivially. The SSOR and the ILU preconditioning amount in this case 
both to the transformation 




(5.2) 



where detM' = detM. To solve a linear system of equations given by M we can now 



apply the Eisenstat trick |35], where we only need to invert triangular matrices. Moreover, 
referring to our criterion mentioned above, we can clearly expect the spectrum of M' to 
be closer to 1 than the spectrum of M. 



^We re-scale the matrix such that the diagonal elements become 1. 
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For the hypercube scalar the situation is more comphcated. To obtain a block struc- 
ture in M, we have to decompose the lattice into 2"^ sub-lattices, which we go through 
successively. This sub-lattice structure is known from staggered fermions. We denote this 
extension of the "red-black" scheme as ^^rainbow preconditioning ^\ Thus we obtain 2'^ unit 
blocks along the diagonal of M, and many more nonzero elements than it is the case for 
the standard action. They arise from the d hopping parameters. We denote the maximal 
magnitude of the elements in L and U as 0(e). We still apply the multiplication (|5.2| ) and 
arrive at 



V V 



i>l 



(5.3) 



We cannot eliminate the off-diagonal elements to such a large extent as it was possible 
for the standard action, because there are more of them from the beginning. On the 
other hand, we can still expect the spectrum of M' to be much closer to 1, typically the 
eigenvalues are 1 — O(e^). Since e is clearly suppressed for the hypercube scalar - compared 
to the £ occurring in the standard action - the suppression of the deviation from the unit 
matrix (0(e) for M, O(e^) for M') becomes even more powerful. In addition we are still 
in agreement with the conditions for the Eisenstat trick. 

To observe this effect, we consider the 2d case for m = 1, where the off-diagonal 
elements are 0, —0.072375 or —0.0213675, as we know from Table ||. On a 6 x 6 lattice one 
obtains the spectrum illustrate in Fig. ^ for M (left), and for M' (right). The hight of the 
lines represents the algebraic multiplicity of the corresponding eigenvalues. We see that 
the eigenvalues of M' are indeed concentrated close to 1, and all of them are < 1, as one 



may suspect from eq. (5^). In particular, the eigenvalue most distant from 1 is 0.72185 
for M, and it moves to 0.94060 for M' . Note also that M does not have an eigenvalue 
1, but in M' there are 18 of them (plus 2 at 0.99537, which are hidden in Fig. 0). The 
"squeezing" effect is indeed stronger for the truncated perfect action than for the standard 
action. 



a. 10 - 
1 

5 



0.7 



.& 10 



0.9 



eigenvalues of M 




eigenvalues of M' 



Figure 7: The eigenvalue distribution of the hypercube scalar matrix M - before precondi- 
tioning - (left), and of the corresponding matrix M' - after preconditioning - (right) on 
a 6 X 6 lattice. 



16 



We conclude that preconditioning is still applicable and efficient for hypercube actions. 
The same holds for hypercube fermions (perfect fermions truncated to a unit hypercube); 
for first experiments in QCD, see Ref. [ |26| . 

6 Non hypercubic lattices 
6.1 Anisotropic lattices 

It is straightforward to write down the perfect free scalar propagator on arbitrary hyper- 
rectangular lattices. We call the lattice spacing in fi direction a^, fi = 1, . . .d. Then the 
blocking from the continuum involves an integral over the hyper-rectangular lattice cells. 
The smearing parameter a is also affected, since it has dimension mass~'^. We have to 
make a choice, for which Id projection we want to obtain ultralocality. If this is the case 
for the mapping on the ly axis, then the generalized perfect propagator reads 

Ganiik) = ^ A(A; + 27r(Va))n2(A; + 27r(//a)) + ^ 

{l/a) := {li/ai,...,ld/ad) 

m) = k,:=^sin^, (6.1) 

where k G B = ] — 7r/ai,7r/ai] x . . . x ] — 7r/ad,iT/ad], and a is still the mass dependent 
smearing parameter from eq. ( ^.17| ). Note that the Poincare invariance of the spectrum 
persists, whatever the lattice spacings. 

Anisotropic lattices are often used in thermodynamics. To simulate at finite temper- 
ature on isotropic lattices, there are typically just Nt = 4: to 8 points available in the 
Euclidean time direction. In order to measure correlators precisely one would like to have 
a larger Nt . The anisotropy permits to invest a larger fraction of the available lattice points 
in the 4 direction (in a hypercubic volume). One introduces a spatial lattice constant a^-, 
and a temporal one at < a^. 

If one then measures masses in all directions, one recovers isotropy up to lattice arti- 
facts, which can be used as a test of improved actions in this context. If the improvement 
is successful, one may use a coarse spatial lattice, but a small at- The dispersion relations 
are the same as on the hypercubic lattice, up to a rescaling of the axes. Hence the strong 
improvement that we observed for the (truncated) perfect scalars persists. 

In improved actions, the presence of "ghosts" (higher branches in the dispersion rela- 
tion) can cause some trouble because they distort the Hermiticity of the transfer matrix. 
In the fully perfect action there are necessarily an infinite number of such ghosts - corre- 
sponding to the sum over I. In Symanzik improved actions there are usually a few ghosts, 
and by some modifications it is possible to get rid of some of them or all of them, if one 
takes in account other disadvantages (for "D234 fermions" this has been discussed in Ref. 
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|28|). The hypercube scalars - as well as all hyper-rectangular scalars - have the virtue 
that they are free of ghosts to start with, cf. Fig. ^. 



6.2 Triangular lattices 



By blocking from the continuum we can also derive efficiently perfect actions on triangular 
lattices or honey comb lattices. Such actions are completely unexplored so far. As an 
illustration we consider the case of a 2d lattice built from regular triangles with lattice 
spacing 1. Here the blocking from the continuum involves integrals over the hexagonal 
lattice cells, which form the dual lattice. f\ In coordinates of the suitable axes - with an 
angle of vr/S - we find that the 11 function introduced in Section 2 is replaced by 



^tria{k) 



3kik2{ki + k2) 



1 ^1,1, ^2 , , ^ ki + k2 

kl cos — + k2 cos [kl + K2) COS 



(6.2) 



We recognize the exchange symmetry of the axes and the special role of the edges of 
the triangular Brillouin zone. Like the hypercubic 11 function, also Utria only contains 
removable singularities, and in the limit ki,k2 — > it becomes 1, which confirms the 
normalization. 

If we only let /c2 ^ 0, we arrive at 



8 



3kl 



Utria{kl,0) 

Remarkably, this is different from n(fci). 



kl kl kl 

1 — cos \ sin — 

2 2 2 



(6.3) 



The blocking constraint 



(6.4) 



— Hx being the hexagon with unit diameter and center x — takes in momentum space the 
form 



kl 



3 

- ^ ip{ki)Iitria{k, 




(6.5) 



Note that I' does not extend over Z^; it only covers the points where I'l — I2 is an inte- 
ger multiple of 3. These points build again a triangular lattice. They are the centers of 
the Brillouin zones, which have a hexagonal form. The lattice momentum k in eq. ( |6.5D 
is in the Brillouin zone around zero, namely the hexagon with the corners (ib47r/3, 0), 
(0, ib47r/3), (47r/3, — 47r/3) and (— 47r/3, 47r/3). Thus the integral over /c, together with the 

®For the perfect action on a honey comb lattice it is the other way round, the cells to be integrated 
over are triangles in two different orientations. 
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sum over I' , just covers the continuum momentum space. Note that this structure agrees 
with that fact that a plane wave along one axis may have momentum An before being 
caught in periodicity. 



Considering that the scalar product is transformed in this coordinate system, we arrive 
at the perfect triangular propagator 

The further steps - evaluation of the inverse propagator in coordinate space, truncation 
etc. - are straightforward, in analogy to Sections 2 to 4. 

7 B-spline blocking 

We return to the hypercubic lattice and we want to explore some new versions of the 
blocking from the continuum. The hope is to find convolution functions, which have 
advantages over the simple prescription described in Section 2, particularly in view of 
locality. 

First we introduce the wavelet or B spline language in d = 1, in a dialect which is 
suitable for our purposes. The zeroth B spline function is just a 5 function, 

/o(x) = 5{x). (7.1) 

Inductively we proceed from order to + 1 as follows: 



1, . , 1 



fN+i{x) = dy My + -) - fN{y - 7t) = / fN{y) dy. (7.2) 



2' 2 



x+l/2 



x-l/2 



This implies for instance 



= { I olherwL ' ^2^^) = { ^ 



|x| l^l < 1 
otherwise ' •'^^'"^ |^ q otherwise 



M^) = { Ul- I^D' I < 1^1 < i ' etc. (7.3) 
otherwise 

We note a few properties: fN{x) = /Ar(-x), /Ar(±(X)) = 0, fN{x) > 0, fN{x)dx = 1, 

fN{x) is maximal in x = and consists of monomial pieces of order A^. Furthermore we 
observe 

E fN{x + j) = 1 (for ah X, iV > 1). (7.4) 

For the periodic blocking from the continuum on a unit lattice, the means that all contin- 
uum points contribute with the same weight to the lattice variables ("democracy among 
the continuum points" ) . 
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For increasing the functions seem flatter and smoother (although is stifl not 
continuously difl'erentiable) . The support ranges from —N/2 to N/2. 

In d dimensions we simply define fN{x) = Y{'li=i fN{x^)- If we use fN{x) as the 
convolution function for blocking matter fields from the continuum, then it is sensible to 
insert 

fN,^l ■■= fN+i{x^) Yl fN{xu) (7.5) 

for the non-compact Abelian gauge field A^. This just provides the connection of nearest 
neighbor lattice matter field variables, i.e. we integrate over all the straight continuum 
connections of corresponding points in adjacent lattice cells. This generalizes the pre- 
scription applied in Ref. [pT| ]. 

The simplest approach in this framework is the decimation RGT (in coordinate space) , 
which defines the lattice variables as 

4>x = j d% /o(x - vMy) , x£Z^. (7.6) 

However, the corresponding expression for the free propagator in momentum space, 

only converges in d = 1. Still the blocking scheme which arises for gauge fields in this 
approach occurs in the literature |^^. On the other hand, for all > 1 the convergence 
of the sum in the propagator is guaranteed, both, for fjy and for Jn^i- 



What has been done so far in this paper, and in Refs. y, 11], is one step forward to the 
use of /i for the matter fields (see eq. ( p. 11 )), and /i^^ for the gauge fields. This already 
provides a promising degree of locality. 

To explore further possibilities, the obvious next step is to consider higher orders in 
this pattern, 

fdS[f[ fN{xu - Vu)] ^{y). (7.8) 

•' u=l 

In momentum space we find inductively 



and the perfect scalar propagator takes the form 



Mk) = U[f) ' (7-9) 



n2^(A; + 27rn 
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decay coefficient 


/l 


/2 


F2 


Cl 


3.805 


2.996 


3.348 




3.876 


3.159 


3.566 


C3/\/3 


3.397 


3.156 


3.757 


C4/2 


3.258 


3.179 


3.333 



Table 3: The competition for locality: we give the coefficients of the exponential decay, 
|;o(i, 0, 0, 0)1 oc exp{—cii}, . . . ,\p{i,i,i,i)\ oc exp{— C4z}, fitted in the interval i = 1...5. 
The decay coefficients Ci are divided by ^/i, so that the level of rotation invariance is 
revealed too. In that respect, the use of F2, and especially f2, is even better than the 
standard blocking function fi (piece-wise constant). 

For increasing N the peak of fN{k) around k = becomes steeper, in agreement with the 
observation that fN{x) turns less localized. Still the Id propagator can be made ultralocal 
for any order N , if we choose the suitable smearing parameter a^- The first few orders 
require at m = 

1 1 P I P k^ ,^ ^ 

an = , ai = — , ao = , as = 1 etc. (7.11) 

" ' 6 ' 4 120 ' 3 40 5040 ^ ' 

The general form is a polynomial in of order — 1. 

The use of these B-spline functions as they stand makes the action less local for N 
increasing beyond 1, so that we cannot find anything better than = 1 (the case we 
had before) within this scheme. The case = 2 is not that bad, but beyond the level of 
locality decreases rapidly. This suggests that a strong overlap of the blocking functions 
for different sites is unfavorable for locality, which was also confirmed by considering other 
types of functions. 

However, we may squeeze the functions again into the interval [—1/2, 1/2] on each axis 
by just rescaling x and correcting the normalization, Fn{x) = N fi^{Nx). Thus we give 
up the property ( |7.4D . Instead of a democratic treatment of the continuum points we do 
something between that and decimation. It turns out that now the locality is in business, 
at least for 

F2{x) = 2(1 - 2\x\) G(l/2 - \x\) . (7.12) 

We summarize these observations in Table |3[ We recommend to focus on A^ = 2, and 
we denote the corresponding blocking function F2 as Eiffel tower function" . It has the 
virtue that it can be related to a relatively simple blocking from a fine to a coarse lattice, 
such that we reproduce the same fixed point. If A^ increases, it becomes more and more 
difficult to establish such relations; the cells, where the block variables are build, have to 
be larger and larger. However, this property is needed if we ultimately want to combine 
the blocking from the continuum with a subsequent multigrid minimization. Moreover, 
the step beyond A^ = 2 is not profitable in view of locality. 
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p(OOOO) 


/9(1000) 


p(llOO) 


p(lllO) 


p(llll) 


4.09671070 


-0.23791704 


-0.056450784 


-0.021127525 


-0.010154672 



Table 4: The couplings p{r) of the truncated perfect scalar at m = 0, constructed from the 
"Eiffel tower" blocking. 

We truncate the FPA obtained from blocking with the Eiffel tower function F2 again 
to a unit hypercube. The couplings of this new hypercube scalar are given in Table 
Its properties regarding the spectrum and thermodynamics are on the same level as the 
hypercube scalar obtained from the standard approach and given in Table ^; all in all the 
hypercube scalar of Section 3 is slightly superior, but in specific respects the Eiffel tower 
hypercube scalar is better. 

For the Wilson-like fermions we observed a qualitatively similar behavior; in that case 
the Eiffel tower blocking looks even a little better than the blocking with the usual step 
function, see Appendix B. 



8 Conclusions and outlook 

We have constructed perfect lattice formulations for free scalar fields of any mass, and 
we optimized their locality. We also explored new aspects of the RGT improvement 
technique, in particular regarding non-hypercubic lattices and blocking functions, which 
are not piece-wise constant. We then truncated the couplings of a perfect scalar to a short 
range - so that the formulation become applicable in simulations - and showed that a 
drastic improvement persists. This was observed from the dispersion relation as well as 
thermodynamic properties. In view of the numerical treatment we showed that a new 
variant of preconditioning is applicable and powerful, which is also true for truncated 
perfect fermions. 

Simulations of truncated perfect scalars have been performed in d = 1 for the anhar- 
monic oscillator |14]. The 2d truncated perfect scalars of Section 3 have been used by W. 



Loinaz in an ongoing numerical study of the critical coupling in the A(j!)^ model, along the 
lines of Ref. |3^. Improved scalar actions play a role in bosonic spin models, in particular 
the non- linear cj- model |2l|, 39 1. A potential field of application in d = 4 is the Higgs 



model (for a very recent application of Symanzik improvement, see Ref. ^^). Finally 
the truncated perfect Laplacian could be applied also for an improved gauge fixing in the 



spirit of Ref. [41|. 
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this paper, a draft of which has been around for more than two years. The first part is 
based on techniques worked out in collaboration with R. Brower, S. Chandrasekharan and 
in particular U.-J. Wiese. 
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A The perfect propagator under a RGT 

We want to demonstrate the perfectness of the propagator G{k) in eq. ( |2.14 ) by expHcitly 
applying a block factor n RGT. The perfect action is supposed to reproduce itself, up to 
a rescaling of the mass. 

From Ref. we can extract the general recursion relation for the propagator under a 
block factor n RGT (the n — > oo limit of which is given in eq. (p.9|)), 



G\k)\ 



1 



m — 9 



a 



(A.l) 



where I G {0, 1,2, , 
m we obtain 

G' {k)\m 



I ^=1 ^?^sin((A;^ + 27r/^)/2n)' 

l}*^. Inserting now the perfect propagator of eq. (2.14) at mass 



n2 4- Ll - l/n2 V 



sin(A;^/2) 



E 



1 



/G^'* (^ + 27rZ + 27rn/)2/n2 + 



nsin((A;,j + 27r/^)/2?i) 

X 



+ 



n 



2sin((A;^ + 27r/"^)/2n) 



sin(A;^/2) 



an[l + 



(fc^j + 27r/^ + 27rn/^)/n nsin((/i;^ + 27r/^)/2n] 



1 



(n^ — l)n' 



2cZ 



En 

I A* 



sin(A;^/2) 



(fc + 2^[r+n/])2 + (nm)2 
G{k)\ 

nm 



sm{{k^ + 27r/^)/2n) 

U{k + 2tt[1 + nl]f 



(A.2) 



This is the desired result, confirming the claims of Section 2. In the last step we have used 
the identity 

sin(a;/2) n2 



n-l 

i:( 

j=0 



sin((x + 2iTj)/2n) 



Amazingly, I could not find this - or a directly related - identity in any table of formulae. 
It can be demonstrated, however, by inserting 



sm{ny/2) 
sin(y/2) 



(n-l)/2 



= ^ exp{iky) 



fc=-(n-l)/2 

(for odd n the sum runs over half integers) and then summing over j. 
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y-x 


pi{y - x) 


\{y - x) 


(0,0,0,0) 





1.83043692 


(1,0,0,0) 


0.158134654 


-0.0577452799 


(1,1,0,0) 


0.0290077129 


-0.0295719810 


(1,1,1,0) 


0.0108088580 


-0.0161376837 


(1,1,1,1) 


0.00476409657 


-0.00889632880 



Table 5: The couplings of the massless hypercube fermion constructed by using the Eiffel 
tower function F2 as RGT blocking function, and by truncating the FPA. 

B Eiffel tower blocking for fermions 

It is straightforward to extend the considerations for non hypercubic lattices and non 
piece-wise constant blocking functions (Sections 6 and 7) to fermions. The most promis- 
ing variant for applications is the Eiffel tower blocking, so we add in this appendix the 
corresponding results for Wilson-type fermions. We denote the lattice Dirac operator as 
D = Pfj,j^ + A, where is odd in p, direction and even otherwise, while the Dirac scalar 
A is entirely even. The prefect free propagator obtained from Eiffel tower blocking - i.e. 



using the blocking function F2 given in eq. ( 7.12 ) - reads 



For massless fermions. Id ultralocality is provided by the RGT smearing parameter ce/ = 
1/2 (this is exactly the same value which is also required for the piece- wise constant 
blocking function /i = Fi |l^]). Here even F7v>2 yields comparable locality of D, but 
for reasons pointed out in Section 7 we prefer to stay with F2. There the locality of the 
FPA is slightly better than for /i: using F2, the decay along the 4d diagonal has the 
exponential coefficients (analogous to C4 in Section 2) 5.037 (for p^J and 5.109 (for A), to 
be compared with 4.931 (for p^) and 5.039 (for A) when we use /i. Therefore truncation 
yields a hypercube fermion, which is at least of the same quality as the one presented 
in Ref. |l^; its dispersion relation does not start with an overshoot, it is excellent up to 
\k\ ~ 2, and also other properties are in business. The new hypercube fermion couplings 
are given in Table |5|. 
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